Hodographic Vortices 

Antonio Moro 
School of Mathematics, Loughborough University, 
Loughborough, Leicestershire, LEll 3TU, UK 
email: a.moro@lboro.ac.uk 



Abstract 

Vortices are screw phase dislocations associated with helicoidal 
wave-fronts. In nonlinear optics, vortices arise as singular solutions 
to the phase-intensity equations of geometric optics. They exist for a 
general class of nonlinear response functions. In this sense, vortices 
possess a universal character. Analysis of geometric optics equations 
on the hodograph plane leads to deformed vortex type solutions that 
are sensitive to the form of the nonhnearity. The case of a Kerr type 
nonlinear response is discussed as a specific example. 

PACS Numbers: 42.15.Dp, 42.65.Hw 
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1 Introduction 

Vortices are fundamental physical objects. They appear in many different 
contexts: from fluid mechanics to nonlinear optics, from superfluidity to 
Bose-Einstein condensation [H El El |4|. In 1974, J.F. Nye and M. Berry 
introduced the concept of phase dislocation obtained by interference of quasi- 
monochromatic wave trains In optics, such phase singularities are called 
optical vortices and, as suggested by Coullet et al. [6], correspond to a state 
analog to vortices in superfluidity. Vortex solitons in Kerr media have been 
numerically predicted and experimentally observed by Swartzlander and Law 
in 1992 [7]. 

The discovery of vortex solitons immediately attracted great interest giv- 
ing birth to a new branch on modern optics referred to as nonlinear singular 
optics [2]. Recent developments, in this field, concern the stabilization of 
optical vortices by propagation in nonlocal nonlinear media [8l[9]. 
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In the linear theory, phase dislocations result from an interference phe- 
nomenon and, consequently, are not observable in the geometric optics limit. 
The situation is completely different in nonlinear regime where the geometric 
optics limit does not exclude the existence of phase dislocations. 

In the present Letter, we are interested in the study of possible new 
singular phase solutions arising in nonlinear geometric optics. Phase dislo- 
cations associated with such solutions have a purely geometric nature and 
their properties are sensitive to the form of the nonlinear optical response. 

We emphasize that the analysis of the long wave limit is also important 
for the study of the full dispersive regime. New singular phase solutions in 
nonlinear geometric optics can be used in the construction of a nontrivial 
ansatz for solving the nonlinear wave equation via, for instance, variational 
or numerical methods. 

The Letter is organized as follows: in section 2, we derive the phase- 
intensity equations for a monochromatic light beam propagating in a generic 
nonlinear medium. We compute the standard vortex (SV) associated with a 
helicoidal wave-front and discuss its universal character. In section 3, a brief 
review of the hodograph method is presented. In section 4, we compute, 
using the hodograph method, a new family of singular phase solutions for a 
Kerr type nonlinear response. We call such solutions hodographic vortices. 
Hodographic vortices turn out to be a geometric deformation of the SVs. 

2 Standard vortex 

Let us consider a stationary paraxial light beam propagating through a 
weakly nonlinear medium of refractive index n^(I) = -|- a^nf(/), where 
/ = |Ep is the intensity of the electric field. We assume to be a 

monotonic increasing function of the intensity, vanishing at / = 0. The con- 
stant no is the value of the refractive index in absence of electric field. The 
weakly nonlinear regime is specified by the maximal value of the intensity 
^max such that the nonlinearity acts as a small perturbation of the back- 
groTind refractive index. Such a perturbation is characterized by the small 
parameter = {v?{Imax) — 'T'o)/'^o associated with the maximal variation 
of the refractive index induced by the electric field. Slow modulations of the 
electric field, propagating along the spatial direction z, are described by the 
following nonlinear Schrodinger (NLS) type equation 




+ konl{\A\^)A = 0, 



(1) 
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where A = A{X, Y, Z) is the envelope of a hnearly polarized electric field 
E = Aexpi{Z/a'^ — ut), = oj/c {c = light speed) and k = kQUQ. 
For our purposes, it is convenient rescaling the equation ([T]) to the dimen- 
sionless form. 

Let us So be the typical spot-size, = ksQ the diffraction length and 
Lnl = ^ / (konKlmax)) the nonlinear length. Introducing the dimensionless 
variables x = X/{^/2so), y = Y/{y/2so), z = Zjil^LniLi), ip = A/y/Imax, 
and the quantity rj = nfLni/riQ, the NLS equation ^ takes the following 
standard form ^ 

iei^, + jV^ + rjM^)tl; = 0, (2) 

where e = Lni/Ld, the subscript denotes the partial differentiation and 
V = {dx,dy). Low dispersion/nonlinear geometric optics limit is obtained 
assuming the diffraction length to be much larger than the nonlinear 
length Lni, i.e. e <C 1. We perform this limit in a standard fashion, looking 
for high oscillating solutions of the form ^p = (pex.p(iS/e). Introducing the 
slow variables u = , v = Sx and w = Sy (note that by definition Vy = Wx) 
it is straightforward to show that, in the limit e — > 0, Eq. ([2]) is equivalent 
to the following dispersionless NLS type equation 

Uz + (uv)^ + {uw)y = (3a) 
Vz + vvx + wvy -rix = (3b) 
■Wz + wwy + vwx — rjy = 0. (3c) 

We refer to the monotonic function 7]{u) as intensity law. We point out 
that, in paraxial approximation, Eq. (|3ap is nothing but the Poynting vector 
conservation law and the equations (|3bH3cl) are equivalent to the eikonal 
equation. 

Let us look for stationary solutions to the system of Eqs. ([3]) such that 
S = z + F{x,y) and u = u{x,y). Note that the functions v = Fx and 
w = Fy are the transverse components of the gradient vector {v, w, 1) which is 
orthogonal to the wavefront (x, y, F{x, y)). In terms of the function F{x, y), 
Eqs. ((3|) are equivalent to the following equations 

F^ + F^ = 2r^{u) (4a) 

U{FxX + Fyy) + UxFx + UyFy = 0. (4b) 

This system of equations is known in fluid dynamics as a model for the 
two-dimensional steady, adiabatic irrotational compressible flow. Function 
F plays the role of the potential velocity and u the density of the fluid. 
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Looking at u as a function of rj (this can always be done due the monotonicity 
of the intensity law), and using the Eq. (l4al) into the Eq. (|4bp . we get the 
quasilinear equation of the form 

AF^^ + BFyy + 2CF^y = 0, (5) 

where A = JF^ + 1, 5 = JF^ + 1, C = JF^Fy and J = d(logu{r]))/dT]. 
The second order equation jS]) is said to be elliptic if its discriminant A = 
AB — = 4Jr/ + 1 is strictly positive. The ellipticity condition A > 
is uniformly (i.e. for any solution) satisfied for few physically relevant 
intensity laws. Important examples are the focussing Kerr-type (u = rp) 
and logarithmic saturable {rj = log(l + u)) nonlinear responses. For this 
reason, in the following, we restrict ourselves to the study of the elliptic case 
only. 

It was observed in Ref. [lO], that, if the function F is harmonic, the 
corresponding wave-front is a harmonic minimal surface. Indeed, provided F 
to satisfy the Laplace equation F^x + Fyy = 0, then F solves the equation ([5]) 
for any J, and consequently for any function u{rf), if and only if it is a solution 
to the minimal surfaces equation 

(1 + F^)F,, + (1 + Fl)Fyy - 2F,FyF,y = 0. (6) 

Moreover, it can be proved that the only harmonic minimal surface of the 
form {x,y, F{x,y)) is the helicoid (x, y, arctan(x/y)) [H]. This singular 
wave-front is the SV associated with screw type phase dislocations [5]. How- 
ever, in nonlinear geometric optics, SVs have a purely geometric origin since 
they do not result from an interference phenomenon. We emphasize their 
'universal' character due to the fact that SVs are solutions to the system of 
Eqs. ((H) for an arbitrary monotonic intensity law. 

In the following, we construct hodographic vortices as a family of vortex- 
type solutions to Eqs. Unlike the SVs, hodographic vortices are not 
universal in the sense specified above. In fact, their geometric structure 
turns out to be depending on the specific form of the nonlinear response. 



3 Hodograph Method 

The system of Eqs. (jlj can be linearized by a hodograph transformation. 
The hodograph method is widely used in fluid dynamics for the study of two 
dimensional compressible flows |12j . 
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Let us introduce the stream function ip via equations 

uv = ipy, uw = -ifx, (7) 

where v = Fx and w = Fy. Using Eqs. ^ into the identity Vy = Wx, one 
gets the following equation for (p: 



u 



{'^xx + ^yy) - Uxfx " Uyify = 0. 



Let us introduce polar coordinates v = pcos9, w = psm9 and suppose that 
the hodograph transformation x = x{p, 0),y = y{p, 0) exists. Expanding the 
total differentials dx{p^ 9), dy{p, 9) and dF{x{p, 9), y{p, 9)), dip{x{p, 9), y{p, 9)), 
we obtain the following set of equations 

cos 9 sin 9 cos 9 sin 9 
-Fp (fp, xe = Fe 



p up ' p up 

sin 9 cos 9 sin 9 cos 9 

Vp = Fp H ifp, yg = Fe H ipg. 

p up p up 

The system of Eqs. ([9]) is over-determined and its compatibility conditions 
Xpd = xgp and ypg = yep lead to the following equations relating F and 

Fe = %p, Fp=p-^(—)^g. (10) 
u op \up J 

The system of Eqs. (fTOl) is also over-determined and it is compatible if and 
only if Fpe = Fep, i.e. 

dp \u^^) ^ dp \up J 

Introducing the variable a = Jqu/p' dp', Eq. pT]l reduces to the following 
generalized Tricomi equation 

+ K {a) •pee = 0, (12) 

where K (a) = —p-§^ i^up)' derivation of Eq. (fT2]) . it is crucial that 

u = u{r]) = u(p^/2) does not depend on the variable 9. 

4 Solution for a Kerr-type nonlinear response 

A focussing Kerr-type medium is specified by the intensity law of the form 
u = co(2?7)'^, where cq and 7 are certain positive constants. In this case, we 
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have a = cop'^'^ / {2^) and u = cop'^^ . Hence, the generalized Tricomi equation 
reduces to the elliptic equation of the form 



fcra H o 



0, 



(13) 



where 7 = (27 + 1) / (47^) . It is a priori not obvious that the exponent of 
a in Eq. (fT3]) does not depend on 7. 

Equation (fT3]) admits the following separable variables solution 



27/ 



(14) 



where a = (1 + \/l + A'yX'^)/2 is a positive real parameter and A is the 
separation of variables constant and plays the role of a geometric deformation 
parameter. Integration of Eqs. (fTO]l gives 



F 



a 



Co \ 2t sin X9 



27 



A 



where a = 27(0 — 1). 

Assuming A^ 7^ 1 and integrating Eqs. 



P 



we obtain 



(15) 



Co 

27 



27 



A{e) /, 



Co 

27 



27 



B{e) p^ 



where 



A{e) 

m 



1 

1 



a sin 9 cos X9 



1 + 27 
27 



A cos 6 sin A^ 



1 



1 + 27 

a cos 9 cos X9 H A sin 9 sin A0 

27 



(16) 



and (3 = 2-f{a - 1) - 1. 

Let us assume that |A| < 1. In this case, it is /3 < and the formula (fTBI) 
implies that p — > 00 as (x, y) —>■ (0, 0). Consequently, the intensity u = cqp^^ 
is divergent at the origin. From the formula (fTSl) it follows that the phase 
S = z + F{x, y) is also divergent at the origin if A 7^ 0. 

We point out that the solution (fTSl) is a deformation of the SV. Indeed, 
we have 



lim F 



arctan {x/y) , 
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On the hodograph plane, the unit vector n = grad S/\gv&d S\ normal to the 
wave front takes the following remarkably simple form 



p 



COS 



e 



p 



sm 



e 



1 



) 



(17) 



Near the origin (i.e. p oo), we have n ~ (cos 0, sin0, 0). Then, the unit 
vector normal to the wave front is completely undetermined. Nevertheless, 
at the point {x,y) = (0,0) the phase (fTSl) is undetermined only for A = 0. 
We call hodographic vortex a solution possessing a point of divergent phase 
where the unit vector normal to the wave-front is completely undetermined. 
From a physical point of view, approaching the origin, the electric field de- 
velops very strong oscillations even for very small As. Moreover, since the 
intensity diverges at the origin, the z— axis is a caustics. Despite the linear 
regime, where caustics and dislocation lines are complementary effects [5], 
in nonlinear geometric optics, vortices generate caustics. Approaching caus- 
tics, geometric optics approximation fails and the wave corrections become 
important. 

Note also that the phase associated with hodographic vortices depends 
on the nonlinearity strength 7. This means that hodographic vortices are 
sensitive to the form of the nonlinearity. 

5 Conclusions 

Hodographic vortices arise in nonlinear geometric optics as deformation of 
SVs. Unlike SVs, the refractive index distribution induced by a hodographic 
vortex is sensitive to the form of the nonlinearity. 

Intriguing features of hodographic vortices stimulate the investigation of 
further significant solutions to the equation (fTSl) . Moreover, the study of 
other physically relevant nonlinear responses will also be of interest. 

We finally stress that screw type vortex solutions studied above arise 
from the analysis of a particular (stationary) reduction of the dispersionless 
NLS type equation ([3|). In general, one can construct an infinite family of 
reductions that are integrable via the hodograph method. Edge and edge- 
screw type phase dislocations are expected to be associated with such non 
stationary reductions. 
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